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Analytical and continuous-time quantum Monte Carlo methods are used to investigate the possi- 
bility of occupation switching and quantum criticality in a model of two quantum impurities coupled 
to two leads. A general discussion of potential occupancy-switching related quantum critical points 
is given, and a detailed analysis is made of a specific model which has been recently discussed. 
For spinless electrons, no phase transition is found. For electrons with spin, a critical value of the 
interaction strength separates a weak coupling regime in which all properties vary smoothly with 
parameters from a strong coupling phase in which occupation numbers vary discontinuously as level 
energies are changed. The discontinuity point is characterized by non-Fermi-liquid behavior. Results 
for self-energies and correlation functions are given. Phase diagrams are presented. 



PACS numbers: 73.21. La, 71.10.Hf, 73.23.-b, 71.27.-fa 



I. INTRODUCTION 

Quantum impurity models (finite interacting sys- 
tems coupled to infinite non-interacting reservoirs) give 
rise to interesting quantum critical behavior including 
the localization-delocalization transition of the macro- 
scopic quantum tunneling problem^i the non-Fermi- 
liquid physics associated with the under- and over- 
screened multichannel Kondo problem and the crit- 
icality related to the two-impurity Kondo modeli^'^'^ 
Recently, a different model apparently exhibiting quan- 
tum criticality was uncovered in the context of exper- 
iments involving a quantum dot inserted into one arm 
of a Aharonov-Bohm interferometer i^i^ As a gate volt- 
age was varied between two Coulomb blockade peaks, the 
transmission phase was found to display a sharp "phase 
lapse" of — TT. It is known that the phase lapse is as- 
sociated with zeroes of the transmission amplitude or 
conductancer^iii and that many-body effects can play an 
important role in producing phase lapses One possible 
mechanism for the vanishing of transmission amplitude 
and thus of "phase lapses" is "population switching" of 
different levels on the quantum dot,i^iiiii^ii^iiIii^ii^ This 
motivates an inquiry into the possibility of obtaining an 
abrupt population switching in models of multilevel dots 
coupled to leads. 

To model the phase lapse system, Golosov and Gefen^S 
introduced a two state quantum dot model, which they 
solved via mean-field theory in the spinless fermion ap- 
proximation. They found population switching which led 
to phase lapse behavior similar to that observed experi- 
mentally but raised the question whether the transition 
survives beyond the mean-field approximation. Medcn 
and collaborators used numerical renormalization group 
methods to show that it in fact does not survive F^^i^ii^ 
In this paper we use analytical and numerically exact 
quantum Monte Carlo (QMC) techniques to go beyond 
mean-field theory and examine the switching transition. 
We find, in agreement with numerical renormalization 



group studiesi^ii^i^ii^ of that in the spinless case studied 
by Ref. [20l . quantum fluctuations destroy the transition, 
whereas in the case of fermions with spin a population- 
switching transition can exist at T = 0. We present 
a brief discussion locating the transition in the general 
landscape of impurity- model quantum phase transitions. 

The rest of this paper is organized as follows. Section 
|TT] presents the model that we study and a general discus- 
sion of the circumstances under which a quantum phase 
transition may take place. Section lllll outlincs an approx- 
imate analytical approach to the problem which follows 
closely Hamann's analysis of the one-impurity Anderson 
model] ^^'^^ Section HVl presents our numerical results and 
Section |V] is a summary and conclusion. 



II. MODEL AND METHODS 
A. Overview 

In this section we present the models to be studied. 
While the initial impetus for the research comes from a 
specific realization of an interferometer involving a quan- 
tum dotiSiiii which implies an impurity model with a 
specific structure?^ it useful to present the results in a 
more general context. 

A general quantum impurity model may be written: 



H = H, 



dot 



lead 
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describing the energetics of electrons in a set of states 
labeled by spin and orbital quantum numbers a, 6, 
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giving the energetics of electrons in a set of infinite leads 
labeled by a momentum (energy) quantum numbers k 
and an index A which denotes lead and spin degrees of 
freedom, and 

H„,i^ = Y,VC^c{,da + h.c. (4) 

kXa 

giving the dot-lead hybridization. The effects of i/ioad + 
i/mix may be encoded in the hybridization functionsi 

F'^^Hc.) = 5:^:i^i^. (5) 

We assume that the lead density of states and hy- 
bridization matrix elements are non-vanishing at the 
Fermi level where = so that 

F-\~iu;n ^ 0) ^ -ir'^'sgniojn) + E^' + O (^^^ (6) 

with Eq a fixed energy characteristic of the conduction 
band band structure and hybridization. We imagine in- 
tegrating out energies on the scale of Eq and higher and 
absorbing the resulting finite renormalizations into the 
parameters of H. We also incorporate the i?"'' into the 
dot parameters e""^ and choose the basis which diagonal- 
izes the level broadening matrix F"^. Refs. [20I I 26I observed 
that transmission experiments which involve electrons in- 
jected and detected in particular combinations of leads 
may have phases which depend crucially on the form of 
the off-diagonal elements when F is written in the basis 
which diagonalizes the input and output currents, but 
these effects are not relevant for the considerations of 
this paper. 

We are interested in possible quantum critical behav- 
ior occurring as parameters in H are varied. In all cases 
known to us, quantum critical behavior may be traced 
back to a crossing of eigenvalues of H as some parameter 
is varied. Coupling to the leads may shift the position 
of the crossing (because e""^ e"'' + E"-^) and in addi- 
tion may either promote the level crossing to a nontrivial 
critical point or convert it to a smooth crossover. For ex- 
ample, a singly occupied dot has a spin 5 = 1/2 so that 
in the presence of a magnetic field one term in H is h- a 
with a the usual Pauli matrices in spin space. If h is held 
parallel to a fixed direction (say z) and the magnitude is 
varied from a positive to a negative value, then a level 
crossing occurs. In the generic case in which only one 
orbital is relevant and the coupling to the leads is anti- 
ferromagnetic the Kondo effect converts the level crossing 
into a smooth crossover. However if more than one or- 
bital is relevant or if the Kondo coupling is ferromagnetic 
a form of multichannel or ferromagnetic Kondo criticality 
may occur One consequence of this criticality is a step 
in the spin polarization as the field is tuned from a posi- 
tive value through zero to a negative value. In a doubly 
occupied multiorbital dot, electronic singlet and triplet 
states are possible, and as a Hund's coupling parameter 



is varied through zero the singlet and triplet states may 
cross in energy. If two channels of conduction electrons 
are present, then coupling to the lattice will in the generic 
case convert the level crossing to a crossover, but in some 
situations the level crossing is shifted to a positive (anti- 
ferromagnetic) value of the Hund's coupling so that if ap- 
propriate symmetry conditions^i^ are imposed the cross- 
ing will become the " Jones- Varma" critical point 

These examples make it clear that quantum criticality 
is associated with the presence of a degeneracy point in 
the Hamiltonian of the isolated level and with a depen- 
dence of the state of the system on the direction from 
which the degeneracy point is approached. Further, one 
sees that a high degree of symmetry must be enforced to 
obtain level crossing or quantum critical behavior, and 
that in particular the possible paths which the system can 
take through parameter space must be constrained. For 
example, in the Kondo example, if the field passes from 
h = hoz to h = —hoz by rotation, then the level crossing 
and any multichannel critical behavior are avoided. 

B. Specific Models 

In this paper we consider two specific models. The 
first, studied by Golosov and Gefen2S is a two- level spin- 
less fermion model. Adopting a pseudospin notation for 
the two orbitals of the impurity model we may write 

H^P;"'^^^ = A-f+L/nina (7) 

with T the triplet of Pauli matrices acting on the orbital 
subspace and A a generalized crystal field splitting. We 
also write for the level broadenings 

r-('o"r-) 

The model is essentially the Anderson impurity model 
with a spin-dependent hybridization. Degeneracy occurs 
in the one electron subspace as the effective crystal field 
splitting is tuned through zero. We will study the model 
assuming that the degeneracy point is approached along 
the direction A || z. For more general directions of A 
our numerical algorithm encounters a severe sign prob- 
lem while the analytic theory becomes notationally much 
more complicated. However, the existence or not of a dis- 
continuity as the system is tuned across the degeneracy 
point should not depend on how the point is approached. 

We also consider the same model, but for electrons 
with spin. A possible H^ot can be written as: 

-^dot" ^ ('^2,tot - "l,tot) + ?7^^ni,tot"-2,tot 

+ C/"ni,tot(ni,tot - l)/2 + C/^^n2,tot(fi2,tot - l)/2 

4- </cxch"5'l • 5*2 "|- t/pair — hop |0,2)(2,0| 

(9) 

which now has a richer level structure which depends on 
the dot occupancy. (Note that fia.tot = J2cr '^aa'^aa-) In 
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the one or three electron sector there are (counting spin 
and orbital degeneracy) four states. In the generic A 7^ 
case the orbital symmetry is lifted while time reversal 
symmetry would protect the spin degeneracy. If only one 
orbital is relevant one obtains the usual Kondo effect, but 
in the special case A = a multichannel Kondo state may 
become possiblei^ 
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FIG. 1: Two-electron energies calculated as functions of level 
splitting parameter A from J/dot with U^^ = 1.5, C/'^^ = 0.5, 
Jcxch = 0.1, Jpair-hop = 0.1 and (7^^ = 0.5 and 1.5. 

In the two electron sector there are the three members 
of a spin-1 triplet and three singlet states. In the occupa- 
tion number basis we may label the singlet states as |0, 2), 
|1,1) and |2, 0). they are generically non-degenerate. 
Fig. [1] shows possible evolutions of the two electron lev- 
els as the orbital splitting A is varied. A positive Hund's 
coupling Joxch is assumed and two values of the intersite 
Coulomb interaction C/^^ arc shown. The three triplet 
states and the |1,1) singlet state have energy indepen- 
dent of A. For large |A| the lowest two electron state is 
one of |0, 2) or |2, 0). As |A| decreases, the energy of this 
favored state increases and level crossings occur. 

If the intersite Coulomb interaction is not too large 
[dotted Hue, double-dotted-dashed line (orange onhne)], 
there is a level crossing either to a triplet state (if 
■/exch < 0, not shown) or to the |1, 1) singlet (Joxcii > 0, 
shown). As noted above, the coupling to the leads would 
lead to a Kondo-quenching of the triplet states. A level 
crossing to this state would in general become a smooth 
crossover, but if particular symmetry conditions are sat- 
isfied a Jones- Varma fixed point would appear. We 
shall present arguments below indicating that the singlet- 
singlet crossing would become a smooth crossover. 

However, if the intersite Coulomb interaction is suffi- 
ciently large (short dashed blue line, solid red line) then 
the first level crossing will be between the |0, 2) and |2, 0)- 
derived states. The degeneracy would be lifted by the 
"pair-hopping" terms, but in the absence of these terms 
a quantum critical point may ensue. 



We observe from this analysis that a level crossing, 
and hence a quantum critical point, is expected only in 
the limit of vanishing pair-hopping. Physically, Jcxcii 
and Jpair-hop arisc from inter-orbital exchange interac- 
tions. These interactions are likely to be very small in 
the quantum dot case envisioned by Ref. |23 where the 
orbitals correspond to spatially separated regions of the 
quanutm dot. 



C. Methods 

We shall be interested in the d electron Green's func- 
tion, defined on the Matsubara axis as: 

C^'ilUr.) = dTe-'^-^Trda{T)dlm, (10) 



G"'' can be expressed in terms of the hybridization func- 
tion and a self-energy resulting from the interaction term 
as: 

[G'"'(ZC^„)]-1 = iuJnSab - - ^^'^'(-iWrO " ^''\iuJn) 

(11) 

and wc shall mainly present results for S. 

We will also present results for the density-density cor- 
relation function W{t) defined as 

W{t) ^ {Tr [ni.tot(r) - n2,tot(T)] [rii,tot(0) - n2,tot(0)]) 

(12) 

To study the models analytically we follow the tech- 
niques of Hamann,— who wrote the single-impurity 
Anderson model^i as a functional integral, decouple the 
interaction with auxiliary fields, identify minima and con- 
sider the action associated with tunneling paths between 
minima. This maps the problem onto a macroscopic 
quantum mechanical model-one dimensional Coulomb 
gasS^for which the key issue is the magnitude of the 
change in scattering phase shift between minima. 

To study the models numerically we use a continuous- 
time QMC technique J ^^1^^ Because we focus on a model 
without complicated terms such as intradot hopping and 
exchange interactions, it is efficient to use the segment 
representation discussed in Ref. [2^. The method encoun- 
ters a serious sign problem if an off-diagonal hybridiza- 
tion function is used, so we employ a basis in which F is 
diagonal. In its current implementation the method re- 
quires that Hdot also be diagonal so we choose e"'' to be 
diagonal in the basis which diagonalizes F. As explained 
above, this does not affect our results. Very recent work 
indicates that the computational cost of dealing with a 
non-diagonal iJdot need not be prohibitive^ so that in 
future work numerical studies of the case [e, F] 7^ may 
be feasible. 

We adopt a semicircular density of states for the con- 
duction bands. 



2*2 



\u\ < 2t 



(13) 
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and choose (fc-indcpcndent) hybridization parameters 
yaX j-gq^ Q-j gQ ^-j^,^^ ^j^g level widths are much less 

than the bandwidth At. The frequency structure of the 
hybridization function affects only non-universal terms 
such as the quantitative location of the critical points. 
We also specialize to the case U^^ ~ U"^"^ to simplify the 
presentation. 

Our calculations were performed on a parallel com- 
puter cluster with 20 dual core 2GHz nodes; a typical 
point requires up to about 15 hours of computer time 
on one CPU. The perturbation orders were typically 
^ 100, but at the lowest temperatures orders up to 
~ 300 were needed. The QMC technique is formulated 
at non-zero temperature: in our work the lowest acces- 
sible temperature is (3t = 800 for calculations of Green's 
functions and self-energies and jSt = 1600 for a single 
point of the correlation function [Figs.HJb) and fTOTb)] . 



III. ANALYTICAL RESULTS 

The spinless fermion model may be viewed as one- 
orbital Anderson model with a spin-dependent hybridiza- 
tion. To study the model analytically we apply the 
methods of Hamann.— i^l Hamann writes the model as 
an imaginary-time path integral, re-expresses the inter- 
action term as ^ [(rii +11.2)^ ^ ("-i ^ "-2)^], decouples 
the ni -|- n2 and ni — n2 interactions with Hubbard- 
Stratonovich fields x and ^ respectively and integrates 
out the fermions. In the situation considered by Hamann, 
the Hubbard-Startonovich field x which couples to ni+n2 
may be absorbed into the chemical potential and its fluc- 
tuations only provide finite renormalizations. In the sit- 
uation of present interest, the field x must be retained. 
After decoupling and integrating out the fermions the 
partition function becomes 



bandwidth limit we obtain 



Z = Zq V£Vx e 



(14) 



with Zq the part of the partition function independent of 
^, x and (in the pseudospin notation of Eq. ([7]) 



S = Tr In [G(7^ + C(t)t^ + ix{T)l] - dr 

Jo 



e(r)2+x(r)2 



U 



(15) 



with the bare level Green function 

Go{iuJn) = 



1 (16) 



and f,r, = i(r"±r22). 

The analysis begins by solving the mean field equations 
dS/d^ = dS/dx = for time independent ^,x. Defining 
ix = r], absorbing a Hartree shift of rj into the chem- 
ical potential and performing the integrals in the wide 
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(17) 



(18) 



with fl the chemical potential shifted by the Hartree term 
from 77. 

For small U there is only one stable solution. For U 
larger than a critical value there is a bifurcation and two 
stable solutions appear. The minimum U occurs at the 
particle-hole symmetric point A = p, ~ where the tran- 
sition is second order. For A or /i 7^ the transition 
shifts to higher U and becomes discontinuous. At the 
particle-hole symmetric point linearizing the equations 
in 77, ^ shows that the critical U is 



[/ph ^ ^^pllp22 



(19) 



while the eigenvector corresponding to the zero eigen- 
value satisfies 



V 



(20) 



To interpret this result we note that (C + 'y) gives the shift 
of level 1 and (^ - 77) the shift of level 2. If T^^ > F" 
then we see from Eq. ((^ that the two mean field solu- 
tions correspond to a larger shift of the broad level and 
a smaller shift of the narrow level. However, what is im- 
portant for the occupation is the level shift relative to 
the level width, and we see that \£, + ri\/T^^ > |C-?7l/r^^ 
so that the relative change in occupation of the narrow 
level is larger. This structure is found in the solution of 
the full nonlinear equations: in the limit of very differ- 
ent level widths the two solutions correspond to states 
in which the population of the narrow level changes sub- 
stantially while that of the broader level shifts much less, 
with the total level occupancy differing in the two stable 
states. An example is given in Fi g. We note, though 
that the important parameter is -^/F^^/F^^ so an extreme 
disparity in level widths is needed to produce a large dif- 
ference in change in occupation. The two solutions in 
general have different energy; for a fixed fi and U there 
is a A (= in the particle-hole symmetric ease) at which 
the energies cross. 

Fluctuations of course alter the mean field predictions. 
In the strong coupling limit we follow Hamann and iden- 
tify the most important fluctuations as "kinks" in which 
the pair ^, 77 tunnels from the neighborhood of one mean 
field solution to the neighborhood of the other. Hamann 
estimates the action associated with these tunnelling 
events from the solution of a singular integral equation. 
The result is an expression for the partition function as a 
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FIG. 2: Electron spectral functions A{ijj) corresponding to 
Hartree-Fock solutions (Eqs ([T7)l . ([T8ll ') for parameters F^^ — 
0.04, F^^ = 0.25, !7 = 0.4 and A = = 0. The units of lo are 
defined by the values of F", F^^ 



sum over kinks occurring at times Tm with a logarithmic 
interaction between them 



Z 
X exp 



dT2 



'To 



A 



i+l 



= 1,3,. 



To 



To 



-1)^+J'ln 



To 



(21) 



with Zo representing the non-tunnelhng contributions to 
the partition function, To a small cutoff time scale of the 
order of the larger of l/F^^'^^, / a tunnelling fugacity of 
order In Utq determined by the bare parameters in the 
problem, A the difference of A from the critical value at 
which the two solutions are degenerate, and, most im- 
portantly, an interaction coefficient K determined by the 
change AJ in scattering phase shifts between the different 
solutions: 



phase and no phase transition is expected, in contradic- 
tion to the results of Ref. [l^. However, if spin is taken 
into consideration the number of channels is doubled and 
K > 2 becomes possible. In the model with spin the 
various interactions require a multiplicity of Hubbard- 
Stratonovich fieldsi^ Solutions to the mean field equa- 
tions involve non-zero values of the decoupling fields cor- 
responding to the dominant interactions, with the other 
decoupling fields being unimportant Assuming that 
U^^ is dominant and neglecting for the moment the pair 
hopping term we write this contribution to the interac- 
tion term as ^ (ni^tot + "-2,tot)^ - ^ ("-i.tot - «2,tot)^, 
decouple the two terms as before and obtain an action 
of the same form as Eq. (|15p except that all of the 
quantities are now to be interpreted as 4 x 4 matrices 
to include the spin degeneracy. The analysis proceeds 
as above except that an extra factor of two multiplies 
the arctan in Eqs. P7)) . as well as the parameter 
A in Eq. (j21|) . The parameter i^, now approximately 

may become 

larger than 2, driving the system into the unscreened 
phase. 

In our work we choose F"/* = 0.04, F^Vi = 0.25, 
implying \Jc = 0.157... for the model with spin (at the 
particle- hole symmetric point). Numerically solving the 
mean field equations we find that at ?7 « 0.32 the change 
Ani| w 0.85 and An2| ~ 0.53 so that the couphng con- 
stant exceeds 2 and we would expect the model to become 
localized. 

As observed above, the pair-hopping term, by mixing 
the two states |0, 2), |2,0) will destroy the phase transi- 
tion, leaving instead a smooth crossover. 

Finally, we note that the |1,1) |2,0) level crossing 
(shown for smaller J7^^ in Fig. [T|), would involve an oc- 
cupancy change of 1/2 electron per spin, in one orbital 
only, so that we would expect the corresponding model 
to be in its screened phase so that no phase transition 
would ensue. 



IV. NUMERICAL RESULTS 
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TT 



(22) 



The second approximate equality applies in the wide 
bandwidth limit where the phase shift gives directly the 
level occupancy. 

It is known22, that for K < 2 the model defined by 
Eq. (j2ip is in its screening phase: tunneling events pro- 
liferate and physical properties vary smoothly as parame- 
ters such as A are varied, whereas for K > 2 the model is 
in the unscreened or localized phase in which tunnelling 
is suppressed and there is a discontinuous change as A is 
tuned across the critical value. 

In the spinless case the maximal change in density in 
either channel for finite C/ is < 1 so the Hamann anal- 
ysis indicates that the model is always in its screening 



A. Spinless fermions 

Fig. [3] shows our numerical results for the c?-occupancy 
as a function of average dot energy in spinless case. We 
took the energy of two levels of the dot to be equal, i.e. 
e^^ = = e. Level one is assumed to be narrower 
(F"/* = 0.04) and level two broader (F^Vi = 0.25). 
The Coulomb interaction U /t = 0.4. The solid and dash- 
dotted lines present the Hartree-Fock results for the level 
1 and 2 occupancies respectively. Within the Hartree- 
Fock approximation there is a first order phase transition 
at e = £* ~ 0.03t - U/2; as e is tuned through e* the 
population of the two levels "switches" abruptly. How- 
ever, the QMC results, indicated by lines with points in 
Fig. [31 shows that at the value of U studied here the 
phase transition is an artifact of the Hartree-Fock ap- 



6 



0.8 



0.6 



0.4 



0.2 




5t=iOO, h. 
5t=100, no 
5t=200, nf 
5t=200, no 
5t=400, nf 
5t=400, no 
HF, nf 
HF, no 



-0.6 



-0.4 



-0.2 



0.2 
(£-hU/2)/t 



0.4 



0.6 



FIG. 3: Main panel: occupancy of level 1 (solid curves) and 
level 2 (dashed curves) of spinless fermion model as function 
of mean level energy e = e^^ = e^"^ shifted by U/2 obtained 
from Hartree Fock (HF, lines without symbols, black on-line) 
and numerically exact QMC calculations (lines with symbols, 
red, blue and magenta on-line). HF is computed at T = 
while QMC is performed at different temperatures as shown 
in the legends. Inset: expanded view of energy dependence 
of dot occupancies in energy regime close to the particle-hole 



symmetric point. Parameters: V^^ /t 
U/t = 0.4. 



0.04, r^7t = 0.25, 



proximation. At half-filling point e = —U/2, the only 
stable solution is (rii) = (n-z) — 1/2, while the weak tem- 
perature dependence (displayed in an expanded scale in 
the inset) indicates that the simulation has accessed the 
low temperature limit of the model. 

As further evidence of the absence of a phase transition 
in the spinless model we present in Fig. [?]the imaginary- 
time density-density correlation function (Eq. (|12p ) com- 
puted at the particle-hole symmetric point for differ- 
ent interaction strengths at various temperatures. Panel 
(a) displays the full imaginary time dependence (note 
W{t) = W{P-t)). In both the intermediate {U/t = 0.4) 
and the strong {U/t — 1.2) interaction case the correla- 
tion function drops as r is increased from t ~ and 
exhibits a minimum at r = /3/2. The value W{t ~ 13/2) 
has a clear temperature dependence which is displayed 
in more detail in Fig.HJb): we see that all curves extrap- 
olate to 0. For <[//<< 0.6 our temperature range 
is sufircient to resolve clearly the behavior expected 
from Fermi liquid theory; for U /t — 0.8, 1.2 the tempera- 
tures numerically accessible to us are not low enough to 
establish a convincing dependence; the extrapolation 
to zero is evident. 

An analytical expression for the Kondo temperature 
for the spinless model has been derived asji^iiii^^i^ 



Tk « 0.2 



^[/(rii + r22) 



■ exp 



TT£{U + e) 

2[/(rii - 



■In 



p22 



(23) 

where the prefactor 0.2 has been set by recovering the 
result of Anderson modeled when taking F^^ = F^^. For 
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FIG. 4: Imaginary time density-density correlation function 
W (Eq. (|12|l 'l computed at putative quantum critical point 
for spinless fermions at parameters indicated, (a) Full imagi- 
nary time dependence of W{t') normalized to its r = value, 
indicating an obvious temperature dependence for both small 
and large f/. (b) Correlation function evaluated at midpoint 
of imaginary time interval, normalized to value at T = O.Oli, 
indicating dependence at low T. Parameters F^^/f — 0.04, 
r^V* = 0.25. 



the 0.4 < U/t < 1.2 studied in this section the Kondo 
temperature varies from Tx/t ~ 0.006 for U/t = 0.4 to 
Tx/t ^ 0.0006 for U/t = 1.2. Our lowest accessable 
temperature. T/t = 1/800, is about twice the U = 1.2 
Tk, explaining why the first stages of the crossover to 
Fermi liquid behavior are evident even in this case. 

Finally, we consider the electron self energy. At the 
particle-hole symmetric point the Matsubara axis self en- 
ergy is purely imaginary and in a Fermi liquid state would 
vanish proportionally to cun as LUn 0. To make the 
analysis more precise we note that Fermi liquid theory 
implies (see e.g. Ref. [35ll36l ) that at low frequency and 
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FIG. 5: Reciprocal of the imaginary part self energy as a func- 
tion of Matsubara frequencies computed for spinless fermions 
at the particle-hole symmetric point for very large interaction 
U/t = 3,4, with r"/t = 0.04, T^V* = 0.25, pt = 100. The 
low frequency downturn is a indication of Kondo bahavior. 



temperature the real-axis self-energy is approximately 



= {I- z- 



2 \n 



^0 



(24) 



with To a scale of the order of the Kondo temperature, 
Z^^ ~ T/Tq, and A a number of the order of unity. This 
form implies that at the particle-hole symmetric point 

E(za;„) = (1 - Z-^)iu;n - i— A^— sgnK) + ... (25) 

Z Iq 

where the ellipsis denotes terms of higher order in tOn. 

Fig. [5] presents the imaginary part of the reciprocal 
of the Matsubara- axis electron self energy S(ia;„) com- 
puted for several temperatures at two very large U val- 
ues. The weak downturn of at a;„/t ~ 0.1 is a signa- 
ture of the crossover to Fermi-liquid behavior. Surpris- 
ingly, a hint of the crossover to the Fermi liquid behavior 
is evident even at temperatures much higher than the 
Kondo temperature. For U/t = 3, the Kondo tempera- 
ture Tji/t = 2 X 10~^, so the beginning of the crossover 
(upturn in l/ImE) can be observed even for T ^ SOOOT/f . 



B. Fermions with spin 

Fig. shows our numerical results for the d- 

occupancy as a function of average dot energy for 
fermions with spin with U^^ = U^^ = 0, the pair hop- 
ping neglected, the level energies and widths as in the 
spinless case, the moderate interaction U^^ = QAt and 
various temperatures. In contrast to the spinless case, a 
clear temperature dependence is evident. The panel (a) 
inset and panel (b) show that near the degeneracy point 
(ni — 712) ~ l/T, so that the dot occupancy exhibits the 




-0.0008 



0.004 



0.008 



FIG. 6: (a) Main panel: occupancy of level 1 (solid curves) 
and level 2 (dashed curves) for fermions with spin as func- 
tion of mean level energy e = e^^ = e^^ shifted by U^^ such 
that the degeneracy point corresponds to e = obtained from 
QMC calculations (red, blue and magenta online), at different 
temperatures as shown in the legends. Inset; expanded view 
of energy dependence of dot occupancies in energy regime 
close to the particle-hole symmetric point, (b) The same data 
plotted as ((n) — 0.5)/(/3t) displaying 1/T temperature depen- 
dence. Parameters: r"/f = 0.04, T^^/t = 0.25, U^'^/t = 0.4. 
Note the x range is different from that of Fig. [S] 



approximately Curie- Weiss behavior expected of a two 
state system with no mixing between the states, in sharp 
contrast to the spinless case. 

At high U^"^, the QMC calculation takes very long time 
to converge, especially at low temperatures. We believe 
that the convergence difficulties arc related to the strong 
suppression of tunnelling between the two states. The 
expense of the calculations means that for the rest of the 
section we restrict ourselves to the particle-hole symmet- 
ric points e = —U^^, and enforce particle- hole symmetry 
by averaging G{t) and G{f3 — r). 

Figure [7] shows the numerical result for the recipro- 
cal of the imaginary part self energy of site 1 on Mat- 
subara axis. The self energy of site 2 is very similar 
and is not shown. Panel (a) is the U^^/t = 0.2 result: 
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FIG. 7: Reciprocal of the imaginary part self energy of site 1 
calculated for fermions with spin as a function of Matsubara 
frequencies at (a) U^'^ /t = 0.2 and (b) U^'^ /t = 0.6 for temper- 
atures indicated, showing qualitative difference in behavior. 
The self energy is the same for spin up and down thus only 
one of them are shown. In panel (a) the dashed lines show 
an extrapolation to zero frequency. Note the x scale are the 
same but the y scale are different. Parameters: Sd ~ —U^^, 
r"/f = 0.04, r^V* = 0.25; 



we see a clear temperature dependence of the low fre- 
quency self-energy. At the relatively high temperature 
Pt = 100 (open circles, black on-line) the extrapolated 
zero frequency limit [ImS](0)/t]~^ is around -50. As tem- 
perature is reduced to (3t — 200 (triangles, red on-line) 
[ImS(0)/t]-i ~ -87. Then at pt = 400 (squares, ma- 
genta on-line) [ImE(0)/i]~^ ~ —160, and at the lowest 
available temperature (31 = 800 (diamonds, blue on-line) 
[lmi;(0)/t]~^ ~ —300. It is clear that the absolute value 
of the intercept at zero frequency increases rapidly as 
temperature is reduced, consistent with a Fermi liquid 
picture in which the ground state is a coherent combina- 
tion of (0,2) and (2,0). 

However, turning to Fig. ^h) we see that at the rel- 
atively larger interaction strength U^'^ /t = 0.6 a clearly 




FIG. 8: Probability of singlet eigenstates as a function of 
interaction strength at particle-hole symmetric point, showing 
dominance of (0, 2) and (2, 0) pair of states for U^^ larger than 



the critical value. 
T = O.Olt. 



Parameters; r"/t = 0.04, T^^ /t = 0.25; 



different behavior occurs. The large difference in y-scales 
between panels (a) and (b) implies that the self energy 
is much larger at U^^ = 0.6. The traces are approx- 
imately linear in frequency and the zero frequency in- 
tercept is small and temperature independent. Fitting 
the low frequency extrapolation gives [lm'E{iu}n)/t]~^ ~ 
—6Aujn — 0.28, indicating that the E ~ iujn expected in 
a Fermi liquid does not occur and that either the ground 
state behavior has been accessed or at least that any 
Kondo scale is far below the lowest measurement tem- 
perature. The non- vanishing intercept is consistent with 
the strong scattering expected if the ground state is an 
incoherent combination of the two valence states. 

It is natural to associate the onset of non-Fcrmi- 
liquid behavior with dominance of the (0,2)/(2,0) pair 
of states, as we have seen that in the strong coupling 
limit tunnelling between these states is suppressed. In 
support of this idea we present in Fig.[8]thc contribution 
of each singlet state to the partition function. In the 
non-interacting case, eigenstate (1,1) has probability of 
1/4 and (2, 0) and (0, 2) each has probability of 1/16. As 
the interaction is increased, the probablity of (1,1) drops 
while that of (2, 0) and (0, 2) increases (and will approach 
1/2 in the infinite interaction limit). Although there can 
be no sharp transition in occupation probabilities in the 
impurity models we consider, the two states (2, 0) and 
(0, 2) become increasingly important as [/^^ is increased 
and are much more important than the (1,1) state in the 
regimes where non-Fcrmi-liquid behavior occurs. 

Figure[5]shows the reciprocal of the imaginary part self 
energy of site one for a finely spaced series of interaction 
at the very low temperature T = t/800. The analysis of 
the spinless fermion case indicated that a weak upturn 
in a plot of this nature indicated a Kondo temperature 
as low as 1/5000 of the measurement temperature. The 
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FIG. 9: Reciprocal of the imaginary part self energy of site 1 as 
a function of Matsubara frequencies for interactions indicated, 
showing change as critical point is approached. Parameters: 
£ = -f/^^ V"/t = 0.04, r^Vi = 0.25; T = t/9>QQ. 

data indicate that while for U^"^ /t = 0.35,0.36 there is 
still some hint of Fermi-liquid behavior, any Kondo scale 
drops extremely rapidly and by U^'^ /t = 0.4 any Kondo 
temperatm-e is lower than 1/(800 x 5000) 3 x IQ-'^t. 

Figure (TU] (a) shows the correlation function W{t) for 
fermions with spin . Comparing to Fig. [3]^ a) we see that 
while the U^"^ /t ~ 0.2 result shown here is rather similar 
to that of U/t = 0.4 shown in Fig.dJ^a), there are funda- 
mental differences between the U^'^/t = 0.6 result shown 
here and that of U/t =1.2 shown in Fig. IHa). The mo- 
tivation of comparing U = 2U^^ is that they produce the 
same Hartree shift, but one should bear in mind that the 
difference in Hamiltonian will lead to a possible difference 
in Tif, although an evaluation of Tk in the spinful case is 
lacking in the literature. The flat range in the middle of 
the U^'^/t = 0.6 and its weak temperature dependence, 
shows strong evidence that fluctuations between the two 
density eigenstates are greatly suppressed. 

The physics can more clearly be seen in Fig. llOf b) 
which shows the temperature dependence of W{t — (3/2) 
at the particle-hole symmetric point. Fermi liquid the- 
ory predicts W{t) ~ l/r^, thus W{(3/2) ~ as T ^ 0, 
while if the spin fluctuation are frozen, W{(i/2) becomes 
a constant at low temperature. The two smallest inter- 
actions shown in Fig. [10] , U^'^ /t = 0, 0.2 (solid line with 
circle, red on-line, long dashed line with square, black on- 
line) clearly reveal the expected behavior. The next 
smallest interaction, U^'^/t = 0.3 (dashed line with up- 
pointing triangle, blue on-line), reveals a downturn sug- 
gestive of a Fermi liquid ground state, although the 
regime is not reached. However the stronger interactions 
U^^/t = 0.4,0.6 (dash-dotted line with down-pointing 
triangle, magenta on-line and dotted line with diamond, 
orange on-line) traces are flat down to the lowest accessi- 
ble temperature T/t = 0.000625. Combining Fig. [7] and 
Fig. (TU] it is clear that a phase transition happens be- 
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FIG. 10: Imaginary time density- density correlation function 
W (Eg I12|l computed at putative quantum critical point for 
fermions with spin at parameters indicated, showing weak T 
and r dependence in non-Fermi-liquid phase and strong T 
and T dependence in Fermi liquid phase, (a) Full imaginary 
time dependence of W{t) normalized to its r = value, (b) 
Correlation function evaluated at midpoint of imaginary time 
interval, normalized to value at T = O.Oli, showing weak 
T dependence in non-Fermi-liquid phase and dependence 
at low T in Fermi liquid phase. Parameters F^^ /t — 0.04, 
r^V* = 0.25. 



twecn U^'^/t = 0.3 and C/^^/t = 0.4 and the system goes 
from a Fermi-liquid like state {U^'^ < U}'^) to a "frozen 
orbital" non- Fermi-liquid state {U^^ > U^^). 

The critical interaction strength U^^ is numerically 
found to be in the range 0.36i < U^^ < 0.39t. This 
is qualitatively consistent with the analytic estimate 
U^^ = 0.32t. A systematic comparison of the Coulomb 
gas prediction of U^^ (based on mean field estimates of 
the phase shifts obtained from the solutions of Eqs. pT|) 
and (HH)) (note that 2J7^^ replaces U in this case) and 
QMC numerical results as a function of p^^/F^^ with 
r^^/t = 0.04 fixed is shown in Fig. [TT] We see a quali- 
tative agreement between the line and the dots; The dif- 
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FIG. 11: The critical interaction strength Uc as a function of 
r^Vr" with r"/t = 0.04 fixed. Red line: result from ana- 
lytical formula. Blue points with error bars: QMC numerical 
resuhs. Parameters: e = -C/", (7" = f/^^ = 0. 
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FIG. 12: Phase diagram of interaction strength U^^ versus 
f/^^ f/" = ?7^^ has been assumed, e" = -f/^^ - U" /2, 
{i = 1,2). Note the offset of zero of x-axis. In the non- Fermi- 
liquid regime a quantum critical point occurs as level energies 
are changed. Parameters: F^V* = 0.04, F^^/t = 0.25. 



ference may come from the approximations we made in 
the mapping to Coulomb gas: We approximated the in- 
tegration path by a hopping between two r-independent 
minima found by minimizing the potential part of the 
action functional, and the hopping path are assumed to 
be linear with width tq, while in reality the hopping path 
is more complicated. 

Fig. [T^ shows a phase diagram indicating the locus of 
criticality in the plane of on-site interaction strength U^^ 
and intersite interaction strength U^'^ at the particle-hole 
symmetric point e" -C/^^ _ U'y2, {i = 1,2). L/" = 
[/^^ has been assumed for simplicity. Above and to the 



left of the two lines the physics is Fermi-liquid like; below 
and to the right a quantum critical point occurs as level 
energies are changed. The regime between the two lines is 
a crossover region in which we cannot determine whether 
it is Fermi liquid within the temperatures studied {T/t > 
1/800). Wc can see that as U^^ increases, the critical U^^ 
needed for a non-Fermi-liquid behaviour increases. 



V. CONCLUSIONS 

In this paper we have studied the possibility of "quan- 
tum criticality" , defined here as a sharp transition in level 
occupancy as the parameters of a multiorbital quantum 
dot are varied. A sharp transition requires a multistabil- 
ity. with more than one locally stable solution, and thus 
in particular requires that a symmetry or physical mecha- 
nism prevents tunnelling between the different potential 
solutions. In the situation of relevance here, the phys- 
ical mechanism is the orthogonality effect arising from 
dissipative coupling to leads. Our work was motivated 
by an interesting proposal arising in the context of a 
two-dot interferometer;^ but we observe that multista- 
bility is of broader interest in the context of potential 
molecular devices. We presented an analysis of the types 
of level-occupancy-related quantum critical points that 
could arise and focussed on the particular situation in- 
troduced by Ref.H^. We used analytical arguments based 
on a mapping to a Coulomb g well as numerical cal- 
culations to show that while the originally studied case 
of spinless fcrmions exhibited only a Fermi liquid behav- 
ior, a model of fermions with spin could in appropriate 
circumstances exhibit a sharp switching. The result for 
spinless electrons is in agreement with previous numer- 
ical renormalization group studies by Refs. [T5lfT6ll2lll2^ . 
(Note that the jumps indicated in Ref. [H arise from a 
multiple-lead multiple-level situation in which the width 
of one level accidentally vanishes.) The key requirements 
arc an intersite Coulomb interaction which is larger than 
the on-site one and the absence of "pair hopping" terms 
in the dot Hamiltonian. This leads to an "orthogonality" 
exponent greater than the critical value of two and hence 
to localization in macroscopic quantum tunneling sence. 
These essential ingedients may be difficult to realize in 
practice, although a strong local electron-phonon cou- 
pling could lead to a polaronic suppression of the on-site 
Coulomb interaction; thus a system that incorporated 
this physics might be an appropriate realization. An in- 
teresting feature of our numerics is that signatures of a 
Kondo effect are visible in the fermion self energy at tem- 
peratures orders of magnitude above the Kondo scale. 

Recently we became aware of a preprint^! reporting 
also that the spinless fermion model does not exhibit 
quantum criticality, and proposing that quantum criti- 
cality could also be realized in a three-lead model (which 
would similarly increase the orthogonality exponent). 
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